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ABSTRACT 


Complete sequences of the plastid gene matK were determined for 62 species of Poaceae from 60 genera, 26 tribes, 
and nine subfamilies to infer phylogenetic relationships. Restio tetraphyllus (Restionaceae) and Joinvillea ascendens 
(Joinvilleaceae) were used as outgroups. Cladistic analysis using PAUP yielded 39 most parsimonious trees with several 
well-supported major lineages. The strict consensus tree shows Streptochaeta and Anomochloa forming the two most 
basal lineages in grasses, followed by Pharus being sister to the remaining species. The other grasses divide into three 
clades: (1) subfamily Bambusoideae (excluding Brachyelytrum) plus Pooideae; (2) Oryzoideae; and (3) subfamilies 
Panicoideae, Arundinoideae, Centothecoideae, and Chloridoideae (termed PACC). Except for Arundinoideae. monophyly 
of each PACC subfamily is generally well supported: however, relationships among subfamilies are unresolved or weakly 
supported. Results obtained using matK sequences are largely consistent with other phylogenies based on molecular 


and structural data, particularly in that relationships among subfamilies remain unclear. 





Interest in the evolution of grasses began early 
in this century with proposed hypotheses based on 
assessment of existing knowledge of the family 
(e.g., Bew, 1929; Hubbard, 1948; Stebbins, 1956, 
1982; Prat, 1960; Clayton, 1981; Tsvelev, 1983). 
Empirical approaches to phylogenetic reconstruc- 
tion of the Poaceae followed those initial hypothe- 
ses, starting with cladistic analyses of morphologi- 
cal and anatomical characters (Baum, 1987; 
Kellogg & Campbell, 1987; Kellogg & Watson, 
1993). Recently, molecular data have provided the 
grounds for phylogenetic hypotheses in grasses at 
the subfamilial and tribal levels. These studies 
were based on information from chloroplast DNA 


(cpDNA) restriction sites and DNA sequencing of 


the rbcL, ndhF, rps4, rpoC2, matK, nuclear ribo- 
somal DNA (nrDNA) 18S and 26S, phytochrome, 
and granule-bound starch synthase genes, as well 
as the noncoding nrDNA Internal Transcribed 
Spacer (ITS) region (Hamby & Zimmer, 1988; Doe- 
bley et al., 1990; Davis & Soreng, 1993; Cummings 
et al., 1994; Hsiao et al., 1994; Nadot et al., 1994; 
Barker et al., 1995, 1999; Clark et al., 1995; Du- 
vall & Morton, 1996; Liang & Hilu, 1996; Mathews 
& Sharrock, 1996; Mason-Gamer et al., 1998; So- 
reng & Davis, 1998; Hsiao et al., 1999). 


Although these studies have refined our under- 


standing of grass evolution at the subfamilial level 
and, to a certain degree, at the tribal level, major 
questions remain to be resolved. Although the basal 
positions of Anomochloeae, Phareae, and Strepto- 
chaeteae have been established, their relative 
placement and taxonomic status are debatable. Un- 
certainties also exist concerning the phylogenetic 
affinities among subfamilies and the taxonomic 
rank of others such as the Oryzoideae. 

In this study, the chloroplast matK gene was cho- 
sen to address these and other questions pertaining 
to higher-level grass systematics. The matK gene is 
~1515 base pairs (bp) in most angiosperms, locat- 
ed within the trnK intron, and functionally may be 
involved in splicing group II introns (Neuhaus & 
Link, 1987; Ems et al., 1995; Hilu & Alice, in 
press a). The effective application of this gene in 
plant systematics (e.g., Johnson & Soltis, 1994, 
1995; Hilu & Liang, 1997; Kron, 1997) and grasses 
(Liang & Hilu, 1996; Hilu & Alice, in press a, b) 
has already been documented. matK is known to 
have relatively high rates of substitution compared 
to other chloroplast genes (see Olmstead & Palmer, 
1994; Johnson & Soltis, 1995). This gene exhibits 
a relatively high proportion of transversions, and 
the 3' region of its open reading frame (ORF) has 





' We thank Nigel Barker, Lynn Clark, Travis Columbus, Jerry Davis, Tarciso Filgueiras, Gary Fleming, Surrey Jacobs, 
David Knepper, A. Nishiwaki, John Randall, Thomas Wieboldt, the Botanical Garden at Bonn, and the Missouri 
Botanical Garden for supplying DNA or plant samples. Seed material for some accessions was kindly provided by the 
U.S. Department of Agriculture, Agriculture Research Service-National Plant Germplasm System. We also thank Ger- 
asimo Borneo, Thomas Borsch. Gerrit Davidse, John MacDougal, and Christoph Neinhuis for their assistance. This 
work was supported by NSF grant &DEB-9634231 and Sigma Xi. 

? Department of Biology. Virginia Polytechnic Institute and State University, Blacksburg. Virginia 24061-0406, U.S.A. 

* Current address: Department of Biology. Western Kentucky University, Bowling Green, Kentucky 42101, U.S.A. 


ANN. Missourt Bor. GARD. 86: 835-851. 1999. 


836 Annals of the 


Missouri Botanical Garden 





MGI, 


|= 


trnK 5' matK trnK 3' 


< < < 
1210R 9R MGI15 


~200 bp 


Figure 1. Diagram of the trnK region including the matK gene. PCR and sequencing primers are indicated with 
arrows. Primer sequences are: MG) = CTACTGCAGAACTAGTCGGATGGAGTAGAT: MG15 = ATCTGGGTTGCTA- 
ACTCAATG; S5-1F = ACCCTGTTCTGACCATATTG; 1210R = GTAGTTGAGAAAGAATCGC; W = TACCCTATCC- 
TATCCAT; 7B = GATTTATCA/GGATTGGGAT; and 9R = TACGAGCTAAAGTTCTAGC, trnK exons and primers 


are not drawn to scale. 


been demonstrated to be quite useful in resolving 
subfamilial, and to a certain degree, tribal relation- 


ships in Poaceae (Liang & Hilu, 1996). 


MATERIALS AND METHODS 
PLANT SAMPLES 


We sequenced the entire matK gene of 62 Po- 
aceae species representing 60 genera, 26 tribes, 
and nine subfamilies (Appendix 1). Subfamilial and 
tribal classification generally follows Clayton and 
Renvoize (1986). Restio tetraphyllus (Restionaceae) 
and Joinvillea ascendens (Joinvilleaceae) were used 
as outgroups because recent studies have demon- 
strated that these two families are closely related 
to grasses (Doyle et al., 1992; Kellogg & Linder, 
1995; Soreng & Davis, 1998, and references there- 
in). 


DNA ISOLATION, POLYMERASE CHAIN REACTION (PCR) 
AMPLIFICATION, AND SEQUENCING 


Leaf tissue was harvested from either green- 
house-grown plants, field-collected plants, or her- 
barium specimens. Total cellular DNA was isolated 
following M’Ribu and Hilu (1996). Because the 
matK gene is part of the trnK intron, we used two 
primers (MG1 or trnK3914 and MG15), located in 
the trnK 5' and 3' exons, respectively, for PCR 
amplification. For sequencing, trnK region PCR 
products were electrophoresed in 0.896 agarose gels 
and DNA fragments of appropriate size excised and 
purified using a QIAquick gel extraction kit (QIA- 
GEN, Inc., Valencia, California). For each acces- 
sion, the entire matK coding region was sequenced, 
utilizing three to six primers (Fig. 1). Sequencing 
reactions were carried out using two different ABI 
Prism" Dye Terminator Cycle Sequencing Kits 
(Perkin Elmer, Norwalk, Connecticut). Most sam- 


ples were electrophoresed in an ABI 373A auto- 
mated DNA sequencer with a stretch gel or in an 
ABI 310 Genetic Analyzer (Applied Biosystems, 
Inc., Foster City. California). Resulting chromato- 
grams were manually edited using Sequence Nav- 
igator 1.0 software (Applied Biosystems Inc., Foster 
City, California). Sequences were deposited in 
GenBank (see Appendix 1). 


SEQUENCE ALIGNMENT AND PHYLOGENETIC ANALYSIS 


Alignment of complete matK sequences was un- 
ambiguous and, thus, done manually. Twelve gaps 
varying in length from 1 to 9 bp were required to 
align sequences (Table 1). Non-random structure in 
the data was tested by using the random trees op- 
tion in PAUP*4.0b2a (Swofford, 1998). The g, val- 
ue for the distribution of tree lengths of 100,000 
random trees was compared using the critical value 
(at a = 0.05) for 500 variable characters and 25 
taxa. Beyond 15 taxa, g, critical values change only 
slightly, allowing them to be used in a conservative 
test with more taxa (Hillis & Huelsenbeck, 1992). 

Phylogenies were generated using Fitch parsi- 
mony as implemented in PAUP, employing heuris- 
tic searches consisting of 1000 replicates of random 
stepwise addition of taxa with MULPARS on and 
tree-bisection-reconnection (TBR) branch swap- 
ping. Gaps were treated as missing data. Sets of 
equally parsimonious trees were summarized by 
strict. consensus. Parsimony-informative gaps are 
mapped onto the strict consensus cladogram (Fig. 
2). Because the transition/transversion ratio (ns/nv) 
in this matK data set is 1.33:1.0, all characters 
were equally weighted. The ns/nv was calculated 
with MacClade 3.01 (Maddison & Maddison, 1992) 
and based on the strict consensus tree. To explore 
the effects of positional weighting, third positions 
of codons were downweighted to one-half that of 
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List of 12 gaps required for alignment of matK sequences of the outgroups Restio tetraphyllus and Joinvillea 


ascendens, and 62 Poaceae species. Gap codes (a-l) are used in the text and Figure 2. Each gap was determined to 


be an insertion (I) or a deletion (D) based on the strict consensus phylogeny (Fig. 2). Position (5' — 3’) is based on 


the overall number of aligned nucleotides (1576). * indicates parsimony-informative indels. 


L———————————————  . . . .. . .. _ CV 
qq  O—,_—PA< —n———— O 


Code Taxa Gap Length (bp) Position 
a Pharus l 6 4-9 
*h All Poaceae D 9 127-135 
*c Leymus. Triticum l 3 154-156 
d Joinvillea, all Poaceae l š 163-165 
e Echinochloa I 3 166-1608 
f Lithachne D 3 283-285 
g Zizania I 6 829-834. 
*h Hyparrhenia, Sorghum l 6 1471-1476 
i Ehrharta D l 1552 
*J All Poaceae except Anomochloa, Streptochaeta D l 1559 
k Restio D 3 1561-1563 
*] PACC I 6 1564-1569 


SEES 


first and second positions. Decay indices (Bremer, 
1988; Donoghue et al., 1992) and bootstrap values 
(Felsenstein, 1985), based on 100 replicates, were 
calculated as measures of support for individual 
clades. Decay analyses were performed with 
AutoDecay (Eriksson & Wikstrom, 1996), using the 
reverse constraint option in PAUP. The data set was 
also analyzed with a Neighbor-Joining (NJ) ap- 
proach utilizing Jukes-Cantor (1969) and Kimura 
2-parameter (1980) distance estimates. 


RESULTS 


MATK LENGTH, GC CONTENT, SEQUENCE DIVERGENCE, 
AND NUCLEOTIDE SITE VARIATION 


The matK ORF in the species examined ranges 
in length from 1521 bp (Ehrharta) to 1548 bp (Hy- 
parrhenia and Sorghum). However, most species 
have an ORF of 1536-1542 bp. Mean guanine + 
cytosine (G + C) content is ~32%. Pairwise di- 
vergence of sequences ranges from 7.3 to 22.3% 
between the outgroups and Poaceae and 0.46 to 
11.9% within Poaceae. Of the overall 1576 aligned 
characters, 836 (53.0%) are variable and 520 
(33.0%) are parsimony informative. In Poaceae 764. 
(48.8%) characters are variable and 473 (30.0%) 
are parsimony informative. The first, second, and 
third positions of codons comprise 30.7, 25.6, and 
43.7% of the variable sites, respectively, and ac- 
count for 28.5, 23.1, and 48.4% of the parsimony- 
informative sites. 


PHYLOGENY OF POACEAE 


Evidence of non-random structure in the Po- 
aceae matK data set is significant (P < 0.01) based 


on the g, value (—0.376). Cladistic analysis with 
gaps coded as missing data yielded 39 equally par- 
simonious trees 2163 steps in length (strict con- 
sensus in Fig. 2). The Consistency Index (CI) and 
Retention Index (RI), excluding uninformative 
characters, are 0.453 and 0.700, respectively. 
Based on the polynomial regression of Sanderson 
and Donoghue (1989), the estimated CI for 64 taxa 
is 0.364; therefore, levels of homoplasy in our data 
set are lower than predicted. Their regression is 
limited to 60 or fewer taxa, and expected ClIs for 
44 and 60 taxa are virtually identical (0.344 and 
0.347, respectively). Thus, the expected CI value 
for 64 taxa according to Sanderson and Donoghue's 
(1989) regression may be inflated. 

The strict consensus tree (Fig. 2) shows Strep- 
tochaeta and Anomochloa to be the two most basal 
lineages, followed by Pharus as sister to the other 
grasses. The remaining species assort into three 
well-supported lineages. One clade includes mem- 
bers of Bambusoideae and Pooideae (including 
Brachyelytrum; Bambusoideae), another contains 
Oryzoideae (including Ehrharta, Ehrharteae), and 
a third clade comprises subfamilies Panicoideae, 
Arundinoideae, Centothecoideae, and Chloridoi- 
deae (PACC). Among the three clades, the PACC 
group has the greatest support with a bootstrap val- 
ue of 100% and a decay index of 13. 

The bambusoid-pooid clade is well supported 
(88% bootstrap and decay of 4) and divided into 
two lineages: one corresponding to the Bambuso- 
ideae and the other to Brachyelytrum + Nardus and 
other Pooideae. Within the Bambusoideae. Sasa 
and Phyllostachys (both Bambuseae) form a strong- 
ly supported subclade (10096 bootstrap) as do Par- 
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iana (Parianeae) and Lithachne + Olyra (both Oly- 
reae). However, Chusquea (Bambuseae) emerges as 
an unresolved lineage. Brachyelytrum (Bambuso- 
ideae, Brachyelytreae; sensu Clayton & Renvoize, 
1986) and Nardus (Pooideae; Nardeae) form a clade 
that is sister to the remaining Pooideae; this clade 
is supported with an 81% bootstrap value and a 
decay index of 3. The other pooids appear in a 
trichotomy including: Stipa + Nassella (Stipeae); 
Melica (Meliceae); and a large, well-supported 
clade comprising Bromeae (Bromus) + Triticeae 
(Hordeum, Leymus, and Triticum) in one subclade 
and Poeae (Poa and Vulpia) + Aveneae (Avena and 
Phalaris) in another subclade. The positions of Bri- 
za (Poeae) and Phleum (Aveneae) are unusual in 
that Briza clusters with the Aveneae and Phleum 
with the Poeae. Subfamily Oryzoideae appears 
monophyletic (100% bootstrap and decay 10) with 
Ehrharta sister to Oryza + Zizania. The Oryzo- 
ideae lie sister to the PACC clade, although boot- 
strap (<50%) and decay (1) support is low. 

In a strongly supported PACC clade, Aristida ap- 
pears basal, but bootstrap and decay support is low 
(Fig. 2). Following this group are four unresolved 
lineages: (1) Arundo, (2) Phragmites + Molinia, (3) 
subfamilies Centothecoideae and Panicoideae, and 
(4) Danthonia, Centropodia, + subfamily Chlori- 
doideae. Monophyly of the Centothecoideae is not 
well supported (74% bootstrap and decay of 1), yet 
support for the Zeugites-Lophatherum-Orthoclada 
subclade is very strong (bootstrap value of 100% 
with a decay index of 15). The Panicoideae do not 
appear monophyletic because Loudetiopsis (Arun- 
dinelleae) appears as an unresolved lineage sepa- 
rate from the Centothecoideae and Panicoideae. 
Aside from Loudetiopsis, the panicoid grasses are 
strongly supported with a bootstrap value of 95% 
and a decay index of 4. The Panicoideae divide into 
two subclades: (1) Paniceae (Digitaria, Panicum, 
and Echinochloa); and (2) the apparently paraphy- 
letic Andropogoneae (Zea, Sorghum, Hyparrhenia, 
and Andropogon), with Tristachya (Arundinelleae) 
sister to the latter three genera. Support for the 
monophyly of, and relationships among, the Cen- 
tothecoideae, Panicoideae, and Loudetiopsis is low. 

The association of Danthonia and Centropodia 
(both Arundinoideae) with subfamily Chloridoideae 
has reasonable support (bootstrap value of 62% and 
a decay index of 2) based on matK data. The Chlor- 
idoideae are very strongly supported (100% boot- 
strap and decay of 10) as a monophyletic lineage. 
Within the chloridoid clade are three well-defined 
subclades, including Uniola, Pappophorum, and 
Eragrostis as sister to the other species, Sporobolus 
+ Zoysia, and a third subclade containing the re- 


maining members of the Chlorideae, Eragrostideae, 
and Orcuttieae. The only well-supported structure 
within this latter group is a clade including Astre- 
bla, Chloris, and Microchloa (97% bootstrap and 
decay 5). 

When third positions of codons are downweight- 
ed (tree not shown) there is a general loss of reso- 
lution among the Arundinoideae taxa at the base of 
the PACC clade, Loudetiopsis emerges unresolved 
and separate from the centothecoid-panicoid clade, 
and the Oryzoideae occupy an intermediate posi- 
tion between the Bambusoideae-Pooideae clade 
and PACC. The NJ analysis based on Jukes-Cantor 
and Kimura 2-parameter distances yielded trees 
identical in topology to each other (Fig. 3) and that 
are largely consistent with the parsimony phylogeny 
(Fig. 2). Most differences are matters of increased 
resolution, such as the relationships among PACC 
subfamilies, and the positions of Chusquea and 
Melica. Topological disagreement does exist con- 
cerning the positions of Orcuttia and Tristachya. In 
the parsimony tree Orcuttia is weakly supported 
(bootstrap <50% and decay 2) as sister to Boute- 
loua; however, in the NJ tree Orcuttia clusters with 
Kengia. The difference in the position of Tristachya 
between the parsimony and NJ analyses is minor, 
involving only a switch with the adjacent Sorghum 
concerning the sister relationship to Andropogon + 
Hyparrhenia. 


INSERTIONS AND DELETIONS (INDELS) IN MATK 


The incorporation of 12 gaps was necessary to 
align the matK sequences of Restio, Joinvillea, and 
Poaceae. From the strict consensus tree (gaps cod- 
ed as missing data), seven insertions and four de- 
letions were determined (Table 1). Another 3-bp 
gap (k) in Restio was identified as a deletion based 
on a broader sampling of matK sequences (Hilu & 
Alice, in press a). Five of these indels are parsi- 
mony-informative (three insertions and two dele- 
tions), of which four are synapomorphic. These syn- 
apomorphic indels include the 9-bp deletion (b) 
uniting the Poaceae, the 1-bp deletion (j) distin- 
guishing Streptochaeta and Anomochloa from all 
other grasses, the 6-bp insertion (l) characterizing 
the PACC clade, and the 3-bp insertion (c) present 
in Leymus and Triticum. The 6-bp insertion (h) 
found in Hyparrhenia and Sorghum is homopla- 
sious. The two indels that are not a multiple of 
three (i and j; both are single base deletions) occur 
near the 3’ end of the matK gene and, thus, do not 
have a major impact on protein composition (Hilu 
& Alice, in press a). Among the remaining indels, 
six (a-f) are located within the first 285 nucleotide 
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sites of the 5’ region, one (g) is found at sites 829— 
834, and three indels (h, k, and 1) occur within 
~100 sites of the 3’ end of the gene. 


DISCUSSION 


Analyses of the matK sequences for the 62 grass 
species demonstrate again the high rate of substi- 
tution in this gene. The 764 (48.8%) variable and 
473 (30.0%) parsimony-informative positions con- 
tribute a considerable number of characters for re- 
solving the phylogeny of the Poaceae. This data set 
is considerably larger than the one used in the ex- 
ploratory matK study of Liang and Hilu (1996) in 
which 17 species were analyzed with 583 nucleo- 
tides, of which 87 were parsimony-informative. The 
overall ns/nv is 1.33 for the whole matK gene. This 
ns/nv ratio is lower than the 1.79 value obtained 
by Liang and Hilu (1996) from the 3' region, but 
is in line with the 1.01 ns/nv ratio calculated by 
Hilu and Liang (1997) for the whole matK gene of 
various plant taxa. The relatively large number of 
transversion mutations detected in the whole matK 
gene of Poaceae appears to be a reflection of dif- 
ferential rates of transversion mutations in different 
sectors of the gene (K. Hilu, unpublished data). 


BASAL LINEAGES IN POACEAE 


The Bambusoideae were traditionally considered 
the most ancestral group of grasses, but the pres- 
ence of derived anatomical, vegetative, and some 
reproductive characters led Soderstrom (1981) to 
state that the Bambusoideae are a specialized 
group. Soderstrom (1981) also stated that among 
Bambusoideae the herbaceous Streptochaeta has 
long been regarded as the most primitive grass. 
Kellogg and Campbell (1987) raised the possibility 
of a most basal position for the *herbaceous bam- 
busoids," while Kellogg and Watson (1993) main- 
tained that the Bambusoideae cannot be both basal 
and monophyletic. 

Using data from ndhF sequences and cpDNA re- 
striction sites alone or in combination with struc- 
tural (anatomy, gross morphology, physiology, and 
chloroplast genome structural mutations) charac- 
ters, Streptochaeta + Anomochloa, and Pharus (tra- 
ditionally considered as herbaceous bamboos) 
emerged as basal lineages with Pharus sister to oth- 
er grasses (Clark et al., 1995; Soreng & Davis, 
1998). Duvall and Morton (1996) also cited support 
for a basal position of Anomochloa based on rbcL 
sequences (Streptochaeta and Pharus were not sam- 
pled). Additionally, the ITS study by Hsiao et al. 
(1999) supports the basal positions of Streptochaeta 
and Pharus (Anomochloa was not sampled). Clark 
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and Judziewicz (1996) placed Anomochloa (An- 
omochloeae) and Streptochaeta (Streptochaeteae) in 
subfamily Anomochlooideae and established sub- 
family Pharoideae to encompass the Phareae. 
Although this study shows Streptochaeta as the 
most basal genus in Poaceae, it does not support 
the monophyly of Anomochloa and Streptochaeta as 
shown by Clark et al. (1995) and in Soreng and 
Davis's (1998) molecular and combined analyses. 
Cladistic analysis of 42 structural characters is 
consistent with our matK results (Soreng & Davis, 
1998). Support for the separation of these two gen- 
era as distinct lineages based on matK sequences 
is reasonable (bootstrap value of 7896 and a decay 
index of 3) and, therefore, these data argue against 
monophyly of the Anomochlooideae sensu Clark 
and Judziewicz (1996). Clark and Judziewicz as- 
serted that it is not easy to find anatomical and 
morphological synapomorphies to define this clade. 
The position of Pharus in the matK phylogeny is 
in agreement with all studies that have included 
this taxon. However, in Soreng and Davis's (1998) 
analysis of structural characters, Pharus was un- 
resolved with Eremitis and these appear as sister to 
other grasses. Pharus also appeared distinct from 
all 215 grass genera (including the Bambusoideae) 
in the numerical study of Hilu and Wright (1982) 


that was based on 85 structural characters. 


DEEP BIFURCATION IN GRASS PHYLOGENY 


When Streptochaeta, Anomochloa, and Pharus 
are excluded from consideration, a split of grasses 
into two lineages is evident, yet the composition of 
the groups varied (Davis & Soreng, 1993; Cum- 
mings et al., 1994; Clark et al., 1995). The phy- 
logeny based on cpDNA restriction sites (Davis & 
Soreng, 1993) depicted two lineages, one corre- 
sponding to the Pooideae and the other including 
all remaining grass taxa. The alliance of the Bam- 
busoideae-Oryzoideae clade with the PACC group 
in the second major lineage was unstable, and the 
overall topology of the tree changed considerably 
after the exclusion of one restriction site (Davis & 
Soreng, 1993, figs. 2, 3). In a subsequent cpDNA 
and structural data study, a bifurcation was not ev- 
ident (Soreng & Davis, 1998). Cummings et al. 
(1994) reported a phylogeny based on the rpoC2 
gene. showing two major grass lineages: (1) Pooi- 
deae + Zea (Panicoideae) and (2) Oryzoideae-Pan- 
icoideae-Arundinoideae-Chloridoideae. The tree 
was rooted with spinach and tobacco, two very dis- 
tant taxa. 

Clark et al. (1995) resolved a bifurcation in the 
family following Streptochaeta + Anomochloa and 


840 


Annals of the 


Missouri Botanical Garden 





J 3 


78 


94 


LE 


100 
28 


Figure 2. 


latifolia, Streptochaeta ang 


88 


indicated with bold, lowercase letters: b 


62 


58 


49 


74 


o 
w 


84 7 


100 

10 
100 
26 


98 


95 114 


3 89 


RESTIO 
JOINVILLEA 


Streptochaeta an 
Streptochaeta spi 


Anomochloa 
Pharus 


Brachyelytrum 


Nardus 
Avena 
Phalaris 
Briza 
Phleum 
Poa 
Vulpia 
Bromus 
Hordeum 
Leymus 
ee m 
rac ium 
i ae 
Nassella 
Stipa 
Lithachne 
Olyra 
Pariana 
Chusquea 
Phyllostachys 
Sasa 
ryza 
Zizania 
Ehrharta 
Aristida ads 
Aristida lat 
Arundo 
Molinia 
Phragmites 
Danthonia 
Centropodia 
Bouteloua 
Orcuttia 
Chloris 
Microchloa 
Astrebla 
Kengia 
Aeluropus 
Monodia 
oysia 
Sporobolus 
ragrostis 
Pappophorum 
Uniola 


Chasmanthium 


Lophatherum 
Orthoclada 
pim ites 
ndropogon 
(ipie i nia 
Sorghum 
Tristachya 
ea 
Digitaria 
Echinochloa 
Panicum 
Loudetiopsis 


oeoproog 


eopriosnquieg 


9 
ƏPƏSpIOZÁIO 


oeapropHO[Q.) — Əeəpiourpuniy | 


əəeəoptooəoul]olu2O 


"V 
e 
B. 
Q 
e. 
e 
A 
[0 


Strict consensus of 39 equally parsimonious trees including 62 Poaceae species and the two outgroups 
Restio tetraphyllus and Joinvillea ascendens (shown in uppercase). Aristida ads 


A. adscensionis. Aristida lat = A. 
S. angustifolia, and Streptochaeta spi = S. spicata. Parsimony-informative indels are 
9-bp deletion, c = 3-bp insertion, h = 6-bp insertion, j = 1-bp deletion, 


and | = 6-bp insertion (see Table 1). For indels, arrows symbolize synapomorphies and bars indicate autapomorphies. 
Numbers above branches are bootstrap values > 50% and numbers below branches are decay indices. 
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Figure 3. Neighbor-Joining phylogram based on Juke-Cantor and Kimura 2-parameter distance models. For abbre- 
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Pharus that corresponded to: (1) Bambusoideae. 
Oryzoideae, and Pooideae; and (2) PACC. This split 
led them to coin the term BOP for the first group, 
although the clade was weakly supported. Mathews 
and Sharrock (1996) reported a similar bifurcation 
based on 174 informative sites from combined phy- 
tochrome gene sequence data. However, the Bam- 
busoideae and Oryzoideae were represented only 
by Bambusa and Oryza, respectively. 

The appearance of subfamilies Panicoideae. 
Arundinoideae, Centothecoideae, and Chloridoi- 
deae in one major group was first demonstrated by 
Hilu and Wright (1982) on the basis of a phenetic 
analysis of morphological-anatomical characters. 
Further support for this grouping was provided by 
the protein and immunological studies of Hilu and 
Esen (1988) and Esen and Hilu (1989). This as- 
semblage, which they named PACC, was also evi- 
dent in the cpDNA restriction site study of Davis 
and Soreng (1993). The monophyly of the group has 
since been substantiated by overwhelming molec- 
ular and structural data (Davis & Soreng, 1993; 
Barker et al., 1995; Clark et al., 1995; Duvall & 
Morton, 1996; Liang & Hilu, 1996; Soreng & Da- 
vis, 1998; Hsiao et al., 1999). Additional support 
for the monophyly of the PACC group is evident 
from the 6-bp insertion at the 3’ end of matK that 
is lacking in other grasses and outside the Poaceae 
(Hilu & Alice, in press a). 

The matK data also show a bifurcation in the 
evolution of the Poaceae. Following the sequential 
divergence of Streptochaeta, Anomochloa, and Pha- 
rus, the Bambusoideae and Pooideae form a clade 
supported by an 88% bootstrap and decay of 4, and 
the Oryzoideae appear sister to PACC. Support for 
the latter relationship is very low (<50% bootstrap 
and decay of 1). Therefore, the BOP clade does not 
gain support from this study, nor is it strongly con- 
tradicted. However, the sister-group relationship 
between Pooideae and PACC is not substantiated 
by these matK data. In contrast with the well-de- 
fined monophyletic PACC lineage, available evi- 
dence is inconclusive concerning the phylogenetic 
affinities among the Bambusoideae, Pooideae, and 
Oryzoideae. 


SYSTEMATICS OF GRASS SUBFAMILIES 

Bambusoideae. The emergence of Streptochae- 
ta, Anomochloa, and Pharus as basal lineages in 
the matK phylogeny provides further evidence for 
the polyphyly of the Bambusoideae as previously 
suspected by Kellogg and Watson (1993) and dem- 
onstrated by other researchers (Barker et al., 1995; 
Clark et al., 1995; Duvall & Morton, 1996; Soreng 
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& Davis, 1998; Hsiao et al., 1999). The clade con- 
taining Bambusoideae (excluding Brachyelytrum) 
shows strong support for the monophyly of the her- 
baceous genera (98% bootstrap and 7 decay). With- 
in this clade the olyroid genera, Olyra and Lith- 
achne, form a well-supported lineage sister to 
Pariana (Parianeae). However, the woody bambu- 
soids (Bambuseae) do not segregate into their re- 
spective subtribes sensu Clayton and Renvoize 
(1986). Sasa and Chusquea of the Arundinariinae 
do not appear monophyletic; instead, Sasa forms a 
strongly supported clade with Phyllostachys (Bam- 
businae), and Chusquea is unresolved. The lack of 
support for the monophyly of these two subtribes is 
also apparent in other studies (Clark et al., 1995; 
Soreng & Davis, 1998; Hsiao et al., 1999). These 
studies also indicate a position for Brachyelytrum 
at or near the base of the Pooideae clade. 

The Pooideae clade (excluding Nar- 
dus) has good support and includes the Stipeae that 
emerge as a basal lineage along with the unresolved 


Pooideae. 


Melica (Meliceae). The Stipeae have been incon- 
sistently treated in the Bambusoideae, Pooideae, 
Arundinoideae, and as a distinct subfamily with af- 
filiation to the Bambusoideae (discussed in Bark- 
worth & Everett, 1987). The basal or near-basal 
position of the Stipeae in the Pooideae is in agree- 
ment with other molecular data (Barker et al., 1995; 
Clark et al., 1995; Mathews & Sharrock, 1996; Cat- 
alán et al., 1997; Soreng & Davis, 1998). 

The remaining pooid genera form two major lin- 
eages with Brachypodium as the sister taxon. The 
Triticeae and Bromeae form one subclade, and the 
Aveneae and Poeae comprise the other. This phy- 
logenetic position for Brachypodium is in agree- 
ment with Catalán et al.’s (1997) results and does 
not differ greatly from that in Soreng and Davis 
(1998) where the genus was sister to the Meliceae 
and basal to most of the Pooideae. Brachypodium 
has been placed in the Triticeae, Bromeae, and 
Brachypodieae (Bor, 1970; Harz, 1980; Hilu & 
Wright, 1982; Clayton & Renvoize, 1986; Macfar- 
lane, 1987). The taxonomic uncertainties regarding 
Brachypodium are due to its floret characteristics 
that are intermediate between the Triticeae and 
Bromeae, possession of smaller chromosomes than 
those found in the Triticeae, and occurrence of base 
chromosome numbers of 7, 9, and 10 (see Hilu & 
Wright, 1982). The matK-based phylogeny herein 
supports the tribal rank of Brachypodium. 

The position of Bromus as sister to the Triticeae 
underscores the phylogenetic affinities between 
these taxa. Bromus is generally placed in its own 
tribe, and Clayton (1978) considered it a link be- 
tween the Poeae and Triticeae. The sister relation- 


Volume 86, Number 4 
1999 


Hilu et al. 
Phylogeny of Poaceae 


843 





ship between the Bromeae and Triticeae is congru- 
ent with the study of Soreng and Davis (1998) and 
Catalán et al. (1997). The three genera of Triticeae 
sampled seem monophyletic, and Leymus and Trit- 
icum appear more closely related to each other than 
either is to Hordeum, as evidenced by a synapo- 
morphic 3-bp insertion (Table 1, Fig. 2). The re- 
maining six pooid genera form two lineages: (1) the 
Poa-Vulpia-Phleum clade, and (2) the Avena-Phal- 
aris-Briza clade. Poa and Vulpia belong to the 
Poeae and Avena and Phalaris to the Aveneae 
(Clayton & Renvoize, 1986). The positions of Briza 
(Poeae) and Phleum (Aveneae) in this analysis are 
anomalous. Yet, other studies that have included 
Briza (Soreng et al., 1990; Hsiao et al., 1999) re- 
solve it within the Aveneae. 

The placement of Melica as an unresolved line- 
age at or near the base of the Pooideae is congruent 
with other molecular phylogenies and supports its 
treatment as a separate tribe. In the NJ tree, Melica 
is sister to the Stipeae (Fig. 3). Brachyelytrum, tra- 
ditionally classified in the Bambusoideae, emerges 
in a clade with Nardus (Pooideae; Nardeae) as sis- 
ter group to the remaining Pooideae. Although this 
relationship may be inconsistent with traditional 
classifications based largely on morphology, it is in 
close agreement with recent phylogenetic studies 
(Clark et al., 1995; Catalán et al., 1997; Soreng & 
Davis, 1998; Hsiao et al., 1999). Each of these 
studies presents alternative relationships for Bra- 
chyelytrum and Nardus, but all indicate that one or 


both genera are basal in an expanded definition of 


Pooideae. The overall morphological and anatomi- 
cal affinity between Nardus and the Pooideae has 
been demonstrated by Hilu and Wright (1982). 
Clayton and Renvoize (1986) asserted that the un- 
usual spikelet of the Nardeae gives no clue to its 
origin, but proposed that it would be better treated 
as an early departure from the pooid line before the 
loss of microhairs. The taxa identified as basal to 
the Pooideae could fit Clayton and Renvoize's no- 
tion of odd genera in an evolutionary transition. A 
number of them have one floret per spikelet, a mix- 
ture of pooid and non-pooid characters, and cur- 
rently appear as relics having different geographic 
distributions. This information may point to an ear- 
ly evolution of a pre-pooid group and considerable 
subsequent diversification and geographic radia- 
tion. 

Oryzoideae. The association of Ehrharta (Ehr- 
harteae) with the Oryzoideae (bootstrap 10096 and 
decay index 10, Fig. 2) provides strong evidence 
for an expanded concept of this subfamily. Histor- 
ically, the taxonomic position of Ehrharta (Ehrhar- 
teae) has been disputed. The taxon has been placed 


within the Bambusoideae (Renvoize, 1986; Watson 
& Dallwitz, 1992), Arundinoideae (Ellis, 1987), 
and Oryzoideae (see Hilu & Wright, 1982). In a 
review of the Ehrharteae, Tateoka (1963) concluded 
that the tribe could be placed in or near the Ory- 
zoideae or near the “arundiform” grasses. An arun- 
dinoid affinity of the Ehrharteae is dismissed be- 
cause of its lack of a 6-bp deletion synapomorphic 
to the PACC clade that includes this subfamily 
(Hilu & Alice, in press a). An oryzoid alliance of 
Ehrharta was evident in the numerical analysis of 
Hilu and Wright (1982). In contrast, Soreng and 
Davis's (1998) cladistic analysis of structural char- 
acters does not support the inclusion of Ehrharta 
in the Oryzoideae. The strong support for the Ehr- 
harta-Oryza-Zizania clade using matK sequences 
clearly demonstrates the phylogenetic affinity of 
Ehrharta to the oryzoid grasses. This phylogenetic 
position for Ehrharta is in agreement with the 
ndhF-based phylogeny of Clark et al. (1995), the 
combined cpDNA restriction site and structural 
data analysis of Soreng and Davis (1998), and the 
ITS-based phylogeny of Hsiao et al. (1999). 

The sister relationship of the oryzoid lineage to 
the PACC clade in this matK parsimony tree (Fig. 
2) is not congruent with studies using other data. 
Those studies have variably placed the Oryzoideae 
in an unresolved trichotomy with the Bambusoideae 
and Pooideae (Clark et al., 1995), sister to the 
Bambusoideae (Barker et al., 1995; Hsiao et al., 
1999). or elsewhere (Duvall & Morton, 1996; Ma- 
thews & Sharrock, 1996; Soreng & Davis, 1998). 
However, this sister position of the oryzoids to 
PACC is not well supported (bootstrap <50% and 
decay index 1), and the bootstrap 50% majority- 
rule tree shows subfamily Oryzoideae sister to the 
bambusoid-pooid clade. Differential weighting of 
the codon positions also places the oryzoids in a 
trichotomy with PACC and bambusoid + pooid lin- 
eages (tree not shown). 

The oryzoid grasses have either been recognized 
as a distinct subfamily or included in the Bambu- 
soideae. Analyses of structural data have been in- 
consistent in terms of taxonomic rank of the ory- 
zoids (Hilu & Wright, 1982; Baum, 1987; Campbell 
& Kellogg, 1987; Kellogg & Watson, 1993; Soreng 
& Davis, 1998). The presence of 10—18 kDa (kil- 
odalton) prolamins endorses the affinities between 
oryzoid and bambusoid taxa (Hilu & Esen, 1988), 
but the low immunological cross-reactivities clearly 
demonstrate a high divergence (Esen & Hilu, 
1989). DNA data have shown the oryzoids as a dis- 
tinct entity (Hamby & Zimmer, 1988; Duvall & 
Morton, 1996; Barker et al., 1995, NJ tree; Clark 
et al., 1995; Soreng & Davis, 1998). In this study, 
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the monophyletic oryzoid clade (including Ehrhar- 
ta) is strongly supported by a 94% bootstrap and a 
decay index of 10. This molecular information thus 
strongly supports the treatment of the Oryzoideae 
as a distinct subfamily, although its phylogenetic 
position is unclear. 

Aristideae. This tribe is represented only by Ar- 
istida in this study. The basal position of the genus 
in the PACC clade is weakly supported (bootstrap 
<50% and decay index 1). The taxonomic position 
of the Aristideae in the Poaceae is disputable be- 
cause of unique anatomical features. The tribe has 
been placed either in the Chloridoideae or Arun- 
dinoideae (reviewed in Hilu & Wright, 1982). Car- 
oline and Jacobs (cf. Jacobs, 1987) found differing 
differentiation for the two Kranz sheaths in Aristida 
species from different ecological habitats, leading 
Jacobs (1987) to believe that Aristida is not dis- 
tantly related to the chloridoid grasses. Based on 
prolamin polypeptide size and immunological sim- 
ilarities, Aristida appeared intermediate between 
the Chloridoideae and Arundinoideae (Hilu & 
Esen, 1990, 1993; Esen & Hilu, 1991). The pro- 
lamin profile of Stipagrostis differs from that of Ar- 
istida, and the former shows low immunological af- 
finities to Aristida, grouping with the Chloridoideae 
(Esen & Hilu, 1991). Sequence data from rbcL 
show Aristida and Stipagrostis to be monophyletic 
and sister to the Chloridoideae (Barker et al., 
1995). The chloridoid affinity of Aristida was also 
apparent in the ndhF-based phylogeny (Clark et al., 
1995). Although the distinctness of Aristida in the 
PACC clade is not in question, the phylogenetic 
position of the Aristideae remains unsettled, espe- 
cially given that the tribe has been represented by 
only the type genus in the majority of recent stud- 
ies. Good representation of the tribe and increased 
resolution are essential before a conclusive assess- 
ment of its taxonomic status and phylogenetic po- 
sition can be determined. The group represents a 
heterogeneous assemblage particularly from ana- 
tomical and physiological perspectives; Stipagrostis 
has a Kranz anatomy that differs from the unique 
Kranz pattern of Aristida; Sartidia lacks Kranz 
anatomy. 

Panicoideae. The Panicoideae do not appear 
monophyletic in our study because the position of 
Loudetiopsis (Arundinelleae) is unresolved in the 
parsimony tree (Fig. 2) and is sister to the Cento- 
thecoideae-Panicoideae in the NJ tree (Fig. 3). 
However, the remaining eight genera sampled do 
form a strongly supported clade divided into two 
lineages: one corresponding to the Paniceae (Dig- 
itaria, Echinochloa, and Panicum); and the other 
representing an apparently paraphyletic Andropo- 
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goneae due to the inclusion of Tristachya (Arundi- 
nelleae). Support for the monophyly of the Pani- 
coideae, excluding Loudetiopsis, is very convincing 
(9596 bootstrap and a decay index of 4). Similar 
results were evident in Clark et al.’s (1995) study 
that showed Danthoniopsis (Arundinelleae) nested 
in a clade comprised of Centothecoideae + Thy- 
sanolaena (Thysanolaeneae; Arundinoideae). This 
clade was sister to the Pancoideae. 

Tristachya (Arundinelleae) is strongly nested 
(10096 bootstrap and 14 decay) within the Andro- 
pogoneae clade of Zea, Sorghum, Andropogon, and 
Hyparrhenia (Fig. 2). The presence of spikelets in 
triads characteristic of the Arundinelleae breaks 
down in Tristachya, where paired spikelets of the 
andropogonoid type are found in some species. Tri- 
stachya also emerged within the Andropogoneae in 
the rbcL and rpoC2 (data set II) studies of Barker 
et al. (1995, 1999). Hsiao et al.'s (1999) ITS phy- 
logeny is the only other study that included an 
Arundinelleae representative. In their analysis, 
Arundinella was sister to the Andropogoneae. Thus, 
it seems that the Arundinelleae are not a mono- 
phyletic tribe, but perhaps distributed among at 
least three lineages. This postulate gains support 
from the molecular study of Mason-Gamer et al. 
(1998) in which the Arundinelleae were not mono- 
phyletic. The Arundinelleae share some spikelet 
features with the Andropogoneae and are thought 
to have given rise to the latter tribe (Clayton, 1981). 
The alliance of Tristachya with the Andropogoneae 
is supported by this matK study but not its ances- 
tral position. 

Centothecoideae. The Centothecoideae were seg- 
regated from the Arundinoideae by Clayton (1978); 
however, their phylogenetic position remains unre- 
solved. Based on the rbcL study of Barker et al. 
(1995), the centothecoid Chasmanthium occurred in 
a clade with Thysanolaena (Arundinoideae; Thysan- 
olaeneae) and was separated from the Panicoideae 
by the arundinoid Gynerium. However, the position 
of Gynerium was described as "equivocal." In the 
ndhF study of Clark et al. (1995), the centothecoid 
genera Zeugites and Chasmanthium appeared in a 
clade containing Danthoniopsis (Panicoideae; Arun- 
dinelleae) and Thysanolaena. In this matK study, the 
Centothecoideae appear monophyletic and related to 
the Panicoideae, although support is low and neither 
Danthoniopsis nor Thysanolaena were sampled. 
Apart from Chasmanthium, the remaining Centothe- 
coideae (Zeugites, Lophatherum, and Orthoclada) are 
very closely related (Fig. 2; 100% bootstrap and de- 
cay 15), raising the question about the position of 
Chasmanthium. The centothecoids cannot be includ- 
ed in the Bambusoideae as treated by Watson and 
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Dallwitz (1992) because they possess the 6-bp in- 
sertion unique to PACC (Hilu & Alice, in press a), 
and based on their phylogenetic position in this and 
other molecular studies. A more comprehensive 
sampling of the Centothecoideae, including poten- 
tially related genera such as Danthoniopsis, Gyner- 
ium, and Thysanolaena, are important prerequisites 
for accurately assessing the monophyly of this group 
and its relationships within PACC. Nevertheless, the 
studies that have included centothecoid representa- 
tives, including this one, demonstrate a greater al- 
liance with the Panicoideae than with either the 
Arundinoideae or Chloridoideae. 

Arundinoideae. The Arundinoideae are known 
to be a taxonomically problematic group. The poly- 
phyletic or paraphyletic nature of the arundinoid 
grasses has been proposed on the basis of morpho- 
anatomical characters (Campbell & Kellogg, 1987) 
and rbcL, rpoC2, and ndhF sequence data (Barker 
et al., 1995, 1999; Clark et al., 1995). The Arun- 
dinoideae sample in Barker et al. (1995, 1999) is 
one of the largest among these studies. They pre- 
sented one of 26 most-parsimonious trees that 
showed the arundinoids split between two major 
clades that are supported by bootstrap values of 
only 33% and 4796. The number of most parsi- 
monious trees and the low bootstrap values make 
it difficult to assess the relationships of the arun- 
dinoid taxa. A more recent study based on the rbcL 
gene (Duvall & Morton, 1996) implied monophyly 
of the arundinoids; however, the study included 
only Arundo and Phragmites. 

Our study does not substantiate a monophyletic 
Arundinoideae. The positions of Arundo and Moli- 
nia * Phragmites are unresolved, whereas Dan- 
thonia and Centropodia are closely related to the 
Chloridoideae. The Phragmites-Molinia clade is 
strongly supported with a 9896 bootstrap and a de- 
cay index of 5 (Fig. 2). The affinity between Molinia 
and Phragmites is apparent in Clayton and Ren- 
voize's (1986: fig. 14) diagram of relationships of 
the Arundineae and is congruent with the results 
of Barker et al. (1995, 1999; used syn. Moliniopsis) 
and the combined analysis of Soreng and Davis 
(1998). 

Chloridoideae. 
doideae, which is often disputed, is strongly sup- 
ported by these matK sequence data (Fig. 2) and a 
more comprehensive matK study (Hilu & Alice, in 
press b). The association of Centropodia and Dan- 
thonia with the chloridoid clade is quite intriguing. 
The sister relationship of Centropodia to the Chlor- 


The monophyly of the Chlori- 


idoideae is in agreement with the rbcL-based phy- 
logeny of Barker et al. (1995), as is the position of 
Danthonia. Centropodia is traditionally placed in 


the Arundineae (Clayton & Renvoize, 1986) or con- 
sidered as a danthonoid (see Barker et al., 1995). 
The genus has a well-developed Kranz anatomy 
(Ellis, 1984) that separates it from the Arundineae 
and allies it with the Chloridoideae. Barker et al. 
(1995) indicated that the lack of haustorial syner- 
gids supports its exclusion from the danthonoid 
grasses. The multinerved glumes and lemmas of 
Centropodia represent traits shared with the Pap- 
pophoreae and Uniolinae, members of a basal 
chloridoid lineage. 

The major tribes Eragrostideae and Chlorideae 
do not appear to be monophyletic. The lack of sup- 
port for the Chlorideae and Eragrostideae as dis- 
tinct lineages was also reflected in the morpholog- 
ical-anatomical study of the subfamily (Van den 
Borre & Watson, 1997) and the matK-based study 
of Hilu and Alice (in press b). The emergence of 
the Pappophoreae, Uniolinae, and Eragrostis in a 
basal clade is congruent with Hilu and Alice (in 
press b). Clayton and Renvoize (1986) placed the 
Uniolinae as a basal group in the Chloridoideae, 
and Van den Borre and Watson (1997) demonstrat- 
ed a near-basal position of the Pappophoreae in 
their phylogeny. 

Another noteworthy group includes Sporobolus 
and Zoysia in a well-supported clade that is con- 
sistent with Soreng and Davis (1998) and Clark et 
al. (1995). The presence of the Chlorideae genera 
Astrebla, Chloris, and Microchloa in a strongly sup- 
ported clade is taxonomically sound. A similar as- 
semblage was also apparent in Van den Borre and 
Watson's (1997) morphological study. A compre- 
hensive systematic study of the Chloridoideae is in 
progress (K. Hilu & L. Alice, unpublished data). 


The matK gene provides sequence information 
sufficient for elucidating evolutionary relationships 
among grass lineages. The results of this study 
identify several well-supported clades that are in 
agreement with other recent molecular studies. 
Most evidence points to Streptochaeta and Anom- 
ochloa as representing the most basal grass line- 
ages. However, the question of whether these two 
genera constitute a monophyletic subfamily or two 
distinct lineages remains unanswered due to con- 
flict among data sets. Pharus has an intriguing po- 
sition. Although Pharus corresponds with Strepto- 
chaeta and Anomochloa by its own distinct lineage, 
the genus is strongly separated from most other 
grasses supporting its subfamilial status as pro- 
posed by Clark and Judziewicz (1996). Yet, Pharus 
shares with other grasses the synapomorphic 1-bp 
deletion that is not found in Streptochaeta and An- 
omochloa (Hilu & Alice, in press a). Subfamilies 
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Bambusoideae, Oryzoideae, and Pooideae can be 
considered monophyletic with individual modifica- 
tion. As others have demonstrated, the Bambuso- 
ideae should exclude Streptochaeta, Anomochloa, 
Pharus, and possibly Brachyelytrum, the last men- 
tioned being allied with Nardus and other Pooideae. 
Oryzoid taxa, including Ehrharta, resolve well from 
the Bambusoideae and are strongly supported as a 
monophyletic unit deserving subfamilial rank. The 
Pooideae are also well supported with their defi- 
nition expanded to include Brachyelytrum. Based 
on existing phylogenetic evidence, there is no con- 
sensus with regard to the relationships among Bam- 
busoideae, Oryzoideae, and Pooideae or to their af- 
finities with the firmly established PACC clade. 

Within PACC, the best supported subfamily is 
the Chloridoideae. The Panicoideae (excluding 
Loudetiopsis) and the Centothecoideae are also tax- 
onomically sound groups, although support is low 
for the relationship between Chasmanthium and 
other centothecoid members. However, the Cento- 
thecoideae may be paraphyletic due to the inclu- 
sion of such genera as Danthoniopsis, Gynerium, 
and Thysanolaena nesting within the centothecoid 
clades in other studies. Finally, the Arundinoideae 
are clearly polyphyletic with some elements at or 
near the base of the PACC clade (Aristida), and 
others (Danthonia and Centropodia) closely related 
to the Chloridoideae. 

Most molecular phylogenies of the Poaceae have 
been inferred from chloroplast genome data and, 
thus, may not be considered truly independent. To 
further resolve the systematic relationships in Po- 
aceae and to test existing hypotheses, robust, nu- 
clear-based phylogenetic analyses are warranted. 
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